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Abstract. 

The understanding and modeling of the structure and evolution of stars is based on statistical 
physics as well as on hydrodynamics. Today, a precise identification and proper description of 
the physical processes at work in stellar interiors are still lacking (one key point being that 
of transport processes) while the comparison of real stars to model predictions, which implies 
conversions from the theoretical space to the observational one, suffers from uncertainties in 
model atmospheres. That results in uncertainties on the prediction of stellar properties needed 
for galactic studies or cosmology (as stellar ages and masses). In the next decade, progress is 
expected from the theoretical, experimental and observational sides. I illustrate some of the 
problems we are faced with when modeling stars and the possible tracks towards their solutions. 
I discuss how future observational ground-based or spatial programs (in particular those dedi- 
cated to micro-arc-second astrometry, asteroseismology and interferometry) will provide precise 
determinations of the stellar parameters and contribute to a better knowledge of stellar interiors 
and atmospheres in a wide range of stellar masses, chemical compositions and evolution stages. 
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1. Stellar internal structure and evolution studies: goals and tools 

Major goals of stellar structure and evolution studies are (i) to characterize and de- 
scribe the physics of matter in the extreme conditions encountered in stars and (ii) to 
determine stellar properties (like age and mass) that trace the history and evolution of 
galaxies and constrain cosmological models. To achieve these goals, we rely on numerical 
stellar models based on input physics that integrate the results of recent theoretical stud- 
ies, numerical simulations and laboratory experiments. The models inputs and outputs 
are chosen and/or validated by comparison with accurate astronomical observations. 

Numerical 2 and 3D hydrodynamical simulations of limited regions of stellar interiors 
and atmospheres are now under reach of computers. They provide valuable constraints 
and data for current standard (ID) stellar models: abundances, convection, rotationally 
induced ins tabilities and mixing, magnetic fields, etc. (see e.g. .Asplund 2005, : Talon, 20071 



Zahnll2007l for reviews). In parallel, the physics of stellar plasmas is studied in the labo- 



ratory with (i) fluid experime nts (study of turbulence in rotating, magnetic fluids, etc., 
see e.g. Richard &: Zahnlll999l ). (ii) particle accelerators (nuclear reaction cross sections, 
etc.) and, (iii) the so-called high energy-density facilities (based on high power lasers or 
z-pinches) which aim at exploring the high temperature and high density regimes found 
in stars, brown dwarfs and giant planets to ge t information on the e quation of state 
(EOS), opacities or thermonuclear reactions (see Remington et al.ll2006l for a review). 

Modern ground-based and spatial telescopes equipped with high quality instrumenta- 
tion are in use or under development (VLT-VLTI, JWST, etc.). They provide very accurate 
data which, after treatment, give access to stellar global parameters (luminosity, radius. 
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mass, effective temperature Toft, gravity log 5, abundances, etc.). On the other hand, seis- 
mic data (such as oscillati on frequencies or aniplitud es) are being obtained in velocity 
from the groun d (see, e.g. 'Bedding & Kieldsen| [2007^ an d in phot ometry by the space 
missions MOST (jWalker et al...2003: ) and CoRoT ijlVIichel et"ai1l2006f). In the next decade . 
valua ble observational data are expected. For instance GAIA ( ESA lioOOt iPerrvman et al 



200lh . to be launched in 2011, will make astrometric measurements, at the micro-arc sec- 



ond level together with photometric and spectroscopic observations of a huge number 
of stars covering the whole range of stellar masses, compositions and evolution stages 
while the Kepler mission, to be launched in 2009, will provide the opportunity to mak e 
asteroseismic observations on a wide range of stars (jChristensen-Dalsgaard et al.ll2007[ ). 

In the following, I discuss the different aspects of stellar modeling, the problems en- 
countered and the perspectives. 



2. Stellar models: input parameters and observational constraints 

2.1. Input physics for stellar models 

Stellar model calculation requires a good description of the physical processes at work. 
Microscopic physics (opacities, EOS, nuclear reaction rates, atomic diffusion) are now 
rather well described which improves the agreement between models and observations. 
However difficulties still arise in the modeling of e.g. (i) cold, dense stars (molecular 
opacities, non- ideal effects in the EOS), (ii) advanced stages of evolution (nuclear reaction 
rates) or (iii) stars to be modeled very accurately like the Sun (see Sect. l3.ip . On the other 
hand, despite important recent progress, macroscopic processes (convection, transport of 
chemicals and angular momentum related to differential rotation an d the role of m agnetic 
fields and internal gravity waves) are not fully understood (see e.g. lTalonll2007 ). 



2.2. Inputs from observations and from model atmospheres 

Stellar models calculation involves inputs and constraints derived from observation. Ob- 
servational data (e.g. magnitudes, colors, spectra, light and velocity curves, astrometric 
data) must be treated to get stellar parameters (luminosity, mass, radius, abundances, 
etc.). Model atmospheres are a crucial step in the analysis, for instance to predict fluxes 
in different bands, synthetic spectra or limb-darkening coefficients. They also provide 
outer boundary conditions for the interior model as well as bolometric corrections and 
color-temperature conversions (to convert the outputs of interior models into quantities 
comparable with observations). Recently, model atmospheres have been improved by the 
bringing-in of better atomic and molecular data and of modern computers and algo- 
rithms. The outputs of 2-3D MHD simulations, including NL TE effects, beg i n to be use d 
to derive abundances, T^ff or log 5 from high quality spectra ( Asplund 2005 : Zahnlboo"?! ). 

However, although the internal errors in the determination of stellar parameters are 
becoming quite sma l l for dwarfs, subgiants and giants of spectral types A to K (see 
Table 1 in Lebretoi] 2005 ). large systematic errors remain for cool and hot stars. For 
instance, the systematic errors on the metallicity of metal poor dwarfs and cool giants 
amount to 0.2-0.3 dex, that is ten times more than the typical internal errors, whil e 
differences in Tog-scales can reach 200 to 400 K ()Gustafssonll2004 lAllende Prietoll2006l ). 

Oscillations have been detected in many stars: solar-like, S Scuti, /3 Cephei, 7 Dor, 
Cepheids, RR Lyrae, SPB, WD, etc. They result from the propagation of acoustic pres- 
sure waves or of gravity waves, depending on the mass, evolution stage, chemical compo- 
sition and excitation mechani sm. Valuabl e constraints can be drawn from the frequencies 
or their combination (see, e.g. Christensen-DalsgaardlQSSi : Roxburgh fc Vorontsov..2003l 
and references therein). 
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3. Modeling calibrators 

The rather few stars for which we have strong or numerous observational constraints 
are used as cahbrators, i.e. they serve to vahdate the models and learn on the physics. 
What is learned from calibrators can then be applied to stars with incomplete or less 
accurate observations. This implies to extrapolate to compositions, masses or evolution 
stages not covered by calibrators. Quantities of astrophysical interest such as age, helium 
content or distance scale can then be derived for very large numbers of stars. 

3.1. Solar modeling: AGS05 mixture and seismology 

The solar photospheric abundances have been re-determined recently on the basis of 3D 
radiative- h ydrodynamical model atmospheres including better atomic data and NLTE 



effects fsee I Aspluncill2005D . The new mixture, referred to as AGS05, gives a global metal- 
licity an d abundances of C, N, O smaller by 30-40% than those given by the GN93 
mixture ( Grevesse fc Noelj 1993 ) derived from ID hydrostatic models. As a result, the 



interior opacity is reduced and the solar model no more satisfies helioseismic constraints 
(convection zone depth, sound speed profile, etc.), the problem being especially acute in 
the region between ~0.4 and ~O.7i?0 (upper part of the radiative zone and shear region 
below the convective zone, i.e. tachocline). Several authors have shown that the present 
uncertainties of the input physics of the solar model, in particular th e opacities or the 



atom ic diffusion velocities, can hardly explain the differences (see e.g. iMontalban et al 



On the other hand, it has been suggested that the n eon abundances could be in 



error by a factor of at least 2, but the problem is still open (jGrevesse et al.ll2007l ) 



3.2. Binary system modeling: the RS Cha and a Centauri binary systems 

The modeling of a binary system consists in reproducing the observed constraints under 
the assumption that the two stars have same age and initial composition. This may allow 
to infer the values of the model unknowns: age, initial helium, and physical parameters 
as the mixing-length parameter for convection or overshooting parameter. 

The binary system RS Cha is an interesting eclipsing, SB2 system whose components 
are A-type oscillating stars in a PMS evolution stage correspondi ng to the onset o f the ea r- 
liest reactions of the CNO cycle. The modeling of the system bv lAlecian et al. I (l2007allbl ) 



at the light of new accurate observations (masses, radii, metallicity) has shown that, 
to get an agreement between models and observations, carbon and nitrogen must be 
depleted with respect to their values in the GN93 mixture. The AGS05 mixture fulfills 
this condition but to assess this result and the values of the system age and helium 
abundance derived from the calibration, it is now necessary to further improve the ob- 
servational data. In particular, the [Fe/H] value should be redetermined using 3D model 
atmospheres and it would be valuable to get individual abundances of major elements 
and better seismic data, the present ones being too coarse to provide useful constraints. 

The binary system a Centauri is the closest and best-known one. The observed global 
parameters are accurate (in particular the accuracy on the intcrferomctric radii is of 1%) 
and seismic observations of both stars have allowed to measure the frequencies of several 
low degree p-modes with an accurac y an 0.3-2.0 fj,Hz. Severa l authors have performed 
a calibration of the system (see, e.g. Miglio fc Montalban 20051 and references therein). 



They conclude that it is difficult to find a set of parameters that satisfies simultaneously 
the global and seismic constraints. As a result, the constraints on the physics of the 
models remain loose and age and initial helium abundance are still poorly determined. 
To progress, it would be interesting to better assess the radii and masses by confirming 
the parallax of the system (different values were obtained from the analysis of Hipparcos 
data) and to further improve and enlarge the seismic data that are still scarce and coarse. 
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3.3. Modeling stars in open clusters 

Members of stellar clusters can also serve as calibrators. They can be studied under the 
assumption that they all have the same age and initial composition but different mass. 
The initial helium abundance and age of a cluster can be derived from the comparison of 
a model isochrone with observations of cluster stars in a color-magnitude diagram (see 
Fig. 3 in Lebreton 2005f) . The surface He abundance Y can be derived from the position 
of the lower main sequence (MS) and the age from the MS turn-off. Such studies require 
accurate observations (parallax, magnitude, Toff, [Fe/H]). The metallicity uncertainty 
affects the Y estimate because of the helium-metallicity degeneracy in the H-R diagram. 
Also, uncertainties on Y and on age come from the bad knowled ge of input physics, a s 
envelope or core convection, rotational mixing or atomic diffusion ( Lebreton et al.l2001 ). 

The observation of binary stars in a cluster can provide additional constraints. For 
instance, the position of the low-mass, non-evolved stars in the mass-luminosity plane 
is related to their helium abundance while physics can be constrained if several binaries 
spanning a large mass range can be observed. Again the accuracy on the parallax and 
metallicity are crucial. A study of the Hyades (the only cluster where individual distances 
have been obtained by Hipparcos and masses measured for a few stars) has shown the 
limitation of the method due to the uncertainties on the metallici ty and input physics a nd 
to the small number of stars with accurate mass determination ("Lebreton e t al.l 2001 ). 

Investigations by e.g. lMazumdar k Antial (|200l[ ): lBasu et al.l (2004): P iau et all (|2005h 
have shown that the seismic analysis of different kinds of oscillators should help to probe 
their inner properties as the outer convection zone depth and helium content or the 
convective core boundary and to estimate their mass and age. For instance, in solar- 
type stars, the higher the helium abundance in the convective envelope, the deeper the 
de pression in t h e adia batic index Fi in the region of second helium ionisation. As shown 
bv iBasu et al. (2004), the helium abundance in the envelope of low-mass stars could 
be derived using the signature of this depression in the p-mode frequencies. This would 
require that low degree p-modes are observed with a frequency accuracy of 0.01% and 
that the mass or radius of the star is known independently. 



4. Deriving astrophysical parameters: the example of stellar ages 

In our Galaxy, the ages of A and F stars are crucial inputs for studying the disc 
while those of old metal poor stars and globular clusters provide valuable constraints for 
cosmology. The uncertainty on age depends on many factors (precision of the observed 
position in the H-R diagram, abundances, knowledge of model input physics as convec- 
tion, rotational mixing, atomic diffusion). The case of globular clusters is discussed by 
Chaboyer (these proceedings). I focus here on the ages of A-F stars. 

4.1. Ages of A and F stars and the size of their mixed cores 

A-F stars have convective cores on the MS and may be fast rotators in the S Scuti 
instability strip. Therefore to model these stars, we are faced with difhculties in describing 
central extra-mixing by core convection overshoot, rotationally induced mixing in the 
radiative zone and we have to estimate the effects of rotation on photometric data. 
These processes modify either the stellar models or the position of the observed star in 
the H-R diagram and in turn affect the age determination. 

The efficiency of mixing in the stellar core determines the quantity of fuel available to 
the star with crucial return on its lifetime. Overshooting of the convective cores produces 
an extra-mixing. In model calculation, this extra-mixing is usually crudely parametrised 
with a coefficient aov, the value of which probably depends on mass, composition and 
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evolutio n stage as shown by empirical calibrations based on binari es and MS width obser- 
vations ( Ribas et al. 2000l : Young et al. 2001 : Cordier et al. 2002). Rotationally induced 
mixing can also bring extra fuel to the stellar engine. Recently, progress has been made in 
the current modeling of the transport of ang ular momentum an d chemicals which results 
fr om differential ro t ation (see the review bv lMathis et al.ll2007l ). 

Goupil &: Taioiil ()2002f ) have calculated models of a typical A-star including either 



overshooting or rotational mixing and have found that these distinct pro cesses cannot be 



discriminated in the H-R diagram where they have similar signatures. I Goupil fc Talon 
also showed that if the A-star pulsates as a i^-Scuti star then the signature of the mixing 
process could be seen in the oscillations frequencies, provided enough modes are observed 
and identified. Today we estimate that an uncertainty on ag e ranging from 13 to 24% 
results from the poor knowledge of the inner mixing processes (jLebreton et al. \\l9M)- We 
expect that in the near future, the improvement on the observed H-R diagram positions 
(in particular on the luminosity with micro-arc-second astrometry) and the availab ility 
of precise seismic data will allow to reduce this uncertainty to 3-5% ( Lebretonll2005^ . 



5. Perspectives in the context of micro-arc-second astrometry 

Today very few calibrators are available to probe the physics of stellar interiors. After 
Hipparcos, ~200 stars have distances accurate to better than 1%. Also, the sample of 
stars with masses and radii measured with an accuracy better than 1% remains small. 
Concerning open clusters, Hipparcos provided precise individual distances only in the 
Hyades while binaries have been analysed in the Hyades and the Pleiades only. There 
is no cluster star where solar-like oscillations have been detected. Concerning A-F stars, 
Hipparcos determined the distances of 10'^ A-F stars with an accuracy better than ~10% 
while seismic data have been obtained for quite many stars, but there are often coarse. 

In the next decade, we expect that observations of stars will increase b oth in num- 
bers and in quality. Missions dedicated to globa l astrometry like GAIA ( ESA 2000l : 



Ferryman et al.ll200lh and SIM (jUnwin et al.ll2007l ) will reach a parallax accuracy better 



than 10 micro-arcseconds. In parallel, the measurements of stellar parameters (magni- 
tudes, temperatures, abundances, masses, radii etc.) are expected to be much improved 
thanks to high resolution spectroscopy, interferometry, etc. High-quality seismic data are 
expected from several missions: CoRoT will reach an accuracy of G.l/iHz on frequencies 
for about 50 targets (mainly solar-like oscillators, f3 Cephei and S Scuti) while Kepler will 
enlarge the seismic sample to thousands of stars with a frequency accuracy of 0.1-0.3/iHz. 

5.1. Expected returns from GAIA 

GAIA astrometry will be complemented by photometry and spectroscopy allowing most 
masses and evolution stages to be unprecedentedly documented. The number of calibra- 
tors of stellar physics will be drastically increased and homogeneous global parameters 
will be provided, e.g. magnitudes, masses and abundances. The parallax of 7 10^(21 10^) 
stars will be measured with an accuracy of at least 0.1%(1%) and the mass of stars in 
17 000 binary systems will be obtained with an accuracy better than 1%. About 120 open 
clusters (up to 1 kpc) will be brought to a level of precision better than the one now 
reached in the Hyades and ^10 binaries per cluster will be observed. Parallax measure- 
ments, accurate to better than 0.5%, will be provided for 5 10^ A stars and 3 10^ F stars. 
Furthermore, while -after Hipparcos- direct distances are available with an accuracy bet- 
ter than 12% for only 11 subdwarfs and 2 subgiants, GAIA will provide (i) precise direct 
distances for very large samples of subdwarfs and for all subgiants up to 3kpc and (ii) 
individual distances with an accuracy better than 10% for stars in 20 globular clusters. 
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We will therefore work on very precise H-R diagrams of very large stellar samples with 
complementary data (as mass, radius, detailed abundances and seismic data) for various 
siibsamplcs of stars. The interpretation of these data in the light of future improvements 
on theoretical, numerical and experimental physics will certainly bring further insight in 
the understanding of the stellar interiors and evolution. 
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